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Abstract. We propose an empirical model to determine the form of energy loss 
of charm quarks due to multiple scatterings in quark gluon plasma by demanding a 
good description of production of D mesons and non-photonic electrons in relativistic 
collision of heavy nuclei at RHIC and LHC energies. Best results are obtained when 
we approximate the momentum loss per collision Apx — apx, where a is a constant 
depending on the centrality and the centre of mass energy. Comparing our results 
with those obtained earlier for drag coefficients estimated using Langevin equation 
for heavy quarks we find that up to half of the energy loss of charm quarks at top 
RHIC energy could be due to collisions while that at LHC energy at 2760 GeV/A the 
collisional energy loss could be about one third of the total. Estimates are obtained 
for azimuthal anisotropy in momentum spectra of heavy mesons, due to this energy 
loss. We further suggest that energy loss of charm quarks may lead to an enhanced 
production of D-mesons and single electrons at low pr in AA collisions. 
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1. Introduction 

The vast amount of data collected at the Relativistic Heavy Ion Collider at Brookhaven 
National Laboratory along with those recently collected from the collision of lead nuclei 
at Large Hadron Collider at CERN have led to the momentous discovery of quark 
gluon plasma (QGP). A large part of the effort in the coming years will be devoted 
to the determination of the precise values of the transport properties of the QGP. In 
this context, the energy loss suffered by quarks of different flavours and gluons as they 
traverse the QGP, undergoing collisions and radiating gluons, is a subject matter of 
considerable topical interest. It is most simply demonstrated by a suppressed production 
of hadrons having large transverse momenta in nucleus-nucleus collisions when compared 
to appropriately scaled productions in pp collisions at the same centre of mass energy 
per nucleon. 

It is often suggested that heavy quarks may lose a smaller amount of energy per 
unit length during their passage through QGP compared to light quarks, due to the 
'dead cone effect' [B12]. The experimental results, however, show similar suppression 
for light and heavy mesons [3]. A recent calculation by Abir et al [4] incorporating the 
generalized distribution of gluons in qc — > qcg and gc — > gcg processes with a proper 
accounting of mass of the heavy quarks, leads to a alE/alx for charm quarks which is 
quite similar to those for light quarks at energies > 10-15 GeV. 

The study of charm quark energy loss provides several unique advantages over 
those of light partons. The charm mesons easily stand out in the multitude of light 
mesons. Most of the charm quarks are produced in initial fusion of quarks and anti- 
quarks (qq — > cc) and gluons (gg — > cc), though a small additional production is 
expected [SI El El El E] from multiple scattering between jets, jets and thermalized 
partons, and thermalized partons. It is not yet clearly established if the charm quarks 
thermalize in the QGP |10j . though it is expected that due to their small numbers their 
impact on the bulk properties of the QGP would be negligible. One also expects that 
due to their large mass, charm quarks will not change their direction as they traverse 
the plasma, though they will slow down, making them excellent probes of azimuthal 
dependence of the conditions of the interacting system. As charm quarks have large 
mass, pQCD remains reasonably valid down to lowest px- There is one additional trait 
which should help us in getting flavour dependence of energy loss of heavy quarks. While 
a u or a d or a s quark or a gluon can fragment into one of many mesons or baryons, a 
charm or a bottom quark would mostly fragment into only a D or a B meson or a charm 
or bottom baryon, respectively. Thus a charm quark after losing energy will appear as 
a D-meson having lower energy. We shall see that this would lead to a characteristic 
enhancement of charm mesons or single electrons at low px- 

The paper is organized as follows. Sec. 2 contains expressions for calculations 
of nuclear modification factor- Raa-, azimuthal anisotropy coefficient- t>2, charm 
production, and our model prescription for energy loss. Sec. 3 contains discussion of 
our results, followed by a summary in Sec. 4. 
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Figure 1. (Colour on-line) Raa and V2 for non photonic electrons at RHIC. Top: 
Momentum loss/per collision cx momentum (left) and oc square- root of momentum 
(right). Bottom: Momentum loss per collision= constant (left) and V2(pt) for single 
electrons. 



2. Formulation 

2.1. Charm Production 

The energy loss of quarks and gluons is most easily seen via the suppressed production 
of hadrons measured using nuclear modification factor, Raa'- 

D , v dN A A/d 2 p T dy . . 

{T A A)da pp /d 2 p T dy 

where A^^^ is the hadron production for the nucleus-nucleus system at a given impact 
parameter, Taa is the corresponding nuclear thickness, and o~ pp is the cross-section for 
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the production of hadrons at the corresponding centre of mass energy /nucleon in pp 
collisions. 

One can now calculate 



da. 



pp 



dyidy 2 d?p T 



2x a x h 



(b). 



2 daij(s,i,u) 



+ fi a \x a ,Q 2 )fr(x b ,Q 



(b), 



do~ij(s, u, t) 



dt 



dt 

/(! + *«) 



(2) 



where px and are the momenta and rapidities of produced charm and anti-charm and 
x a and Xb are the fractions of the momenta carried by the partons from their interacting 
parent hadrons. These are given by 



x n 



7? 



+ e 



y-2 ' 



M T 



-yi 



(3) 



where Mt is the transverse mass, ymg+p^, of the produced heavy quark. The 
subscripts i and j denote the interacting partons, and fi/j are the partonic distribution 
functions for the nucleons. The fundamental processes included for LO calculations are: 



9 + 9 -> c + c 

q + q -)■ c + c . (4) 

We recall that the above LO pQCD expression reproduces the NLO results [TTj 
when supplemented with a f^-factor w 2 (see Ref. |12j). 

We have used T^a= 225 fm -2 for 0-10% centrality for Au+Au collisions at RHIC, 
as calculated from Glauber formalism. For Pb+Pb collisions at LHC, T^a =195 fm -2 
for 0-20% centrality has been used. We use CTEQ5M structure function along with 
EKS98 [13J shadowing function. The factorization, renormalization, and fragmentation 
scales are chosen as J tuq + p\ and the charm quark mass has been taken as 1.5 GeV. 

2.2. Energy Loss 

We propose an empirical model for the energy loss for charm quarks which is inspired 
by a multiple scattering model used earlier by [H] supplemented with considerations 
of Baier et al [15] for partonic energy loss [16]. A straight-forward empirical 
implementation of a similar energy loss procedure for light parton has been found to 
provide a satisfactory dependence of energy loss on the centralitiy of the collision |17j . 
We recall once again that the nuclear modification of heavy meson production is similar 
to those for light mesons. 

We perform a Monte Carlo implementation of our model calculations and estimate 
the momentum loss of charm quarks and nuclear modification of D-meson and single 
electron production. We assume that the energy loss of heavy quarks proceeds via 
multiple collisions and that the momentum loss per collision is given by, (see for example 
Ref. PS]) 

(Ap)* = a ( Pl f , (5) 
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so that one can write 

dp Ap 

T x = -- (6) 

where a and (3 are parameters to be determined and A is the mean free path of the 
charm quark, taken as 1 fm, in these initial studies. In what follows, we shall consider 
charm quarks at central rapidities and thus p = pt- Thus the momentum of the charm 
quark after n collisions will be given by 

Pn+1 =Pn~ (Ap) n (7) 

The charm quark can continue to lose energy in collisions as long as the resulting 
momentum remains positive. We estimate the probability for the charm quark to have 
n collisions, while covering the path length L from a Poisson distribution 

P(n,L) = ^e-^\ (8) 
n! 

Taking a value for the coefficient a and the exponent /3, we estimate the largest 
number of collisions- N, which the charm quark having momentum p? can undergo. 
Next we sample the number of collisions n, which the charm undergoes from the 
distribution 

TV 

p{n) = P(n,L)/ y £P(n,L) (9) 

n=l 

to get the final momentum of the charm quark. 

Finally we fragment the charm quarks into D-mesons. Thus we have, 

E % = Eq ^g¥^ ® D{Q Hq) ® f{Hq ~* e) ' (10) 

where the fragmentation of the heavy quark Q into the heavy-meson Hq is described 
by the fragmentation function D. We have assumed that the shape of D(z), where 
z = Pd/Pc, is identical for all the D-mesons, [TH] and 

z[l - 1/z - e p /(l - z)} 2 
where e p is the Peterson parameter and 



o 



dzD(z) = 1 . (12) 



We have kept it fixed at e p =0.13. 

Here F(Hq — > e) denotes semileptonic decay of D-mesons and the electron 
distribution is taken from Ref. [20J. 

2.3. Azimuthal Anisotropy 

Non-central collisions of identical nuclei will lead to an oval overlap zone, whose length 
in and out of the reaction plain would be different. Thus, charm quarks traversing 
the QGP in and out of the plain will cover different path lengths and lose differing 
amount of energy. This would lead to an azimuthal dependence in the distribution of 
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resulting charm mesons, whose azimuthal anisotropy could be measured in terms of the 
i>2 coefficient defined by 

/ f-JZ cos (20) 

V2(PT) = ] T 1 d N [ ( 13 ) 

We have approximated the colliding nuclei as having a uniform density with radius R in 
these calculations and obtained average path-length for the charm quarks along a given 
using Eq. 9 of Ref . [T7] . 
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Figure 2. (Colour on-line) Raa and t>2 for D mesons at LHC. Top: Momentum loss 
per collision cx momentum (left) and oc square root of momentum (right). Bottom: 
Momentum loss per collision^ constant (left) and V2(pt) for D-mesons (right). 
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Figure 3. (Colour on-line) Left: Raa for single electrons for collisions of gold nuclei 
at RHIC and lead nuclei at LHC at b — 4.5 fm. Inset shows enhanced production of 
single electrons having low pt- Right: Raa for D- mesons for the same collisions at 
RHIC and LHC and the inset shows the enhanced production of D-mesons having low 
Pt- The best fit values of a from Figs. [TJ and [2] and (3=1 are used for these calculations. 



3. Results and Discussion 

Let us now discuss our results for nuclear modification factor and azimuthal anisotropy. 
We consider three values for the exponent f}\ 0, 0.5, 1.0 appearing in Eq. [5J inspired by 
the three energy loss mechanisms, namely those applicable in the so-called Bethe-Heitler 
regime, LPM regime, and complete coherence regimes considered by Baier et al. [151 E] 
which lead to energy loss per unit length as proportional to energy, square-root of the 
energy, and independent of the energy for light partons. Kampfer et al [2T] had earlier 
used this approach to study the effect of charm quark energy loss on the correlated 
charm decay. 

Next we vary a to get a description of the Raa for single electrons at RHIC (Fig. 
and for D-mesons at LHC (Fig. [2]). 

In Fig. [T] we show our results for nuclear modification factor Raa and azimuthal 
anisotropy for non photonic electrons at top RHIC energy y/s = 200 A GeV [22J. 
Comparing the results of Figs. [U we see that the model assuming momentum loss 
per collision as proportional to the momentum closely follows the shape of the 
experimentally determined Raa for single electrons almost over the entire range of pt 
under consideration. We add that our theoretical calculations have not included the 
b — > e contribution which is contained in the experimental results which can modify the 
Raa for larger px by up to 10% as the produced b quarks are much less in number and 
also lose much smaller energy |23j. The scenario, where Ap oc */p, is only moderately 
successful in describing the data over a limited pr range of 2-4 GeV/c (Fig. [p. While 
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the assumption of a constant momentum loss per collision may bracket the Raa over 
the very limited range of 2-3 GeV/c, it does not follow the shape of the px dependence 
(Fig. II). 

The best values of the a determined from the results in Fig. [1] are used to estimate 
V2 for the single electrons. We see that our calculations provide a reasonable description 
of V2(pt) f° r Pt > 2 GeV/c and overestimate the results for lower p?. A relaxation of 
our assumption of a static medium at a constant temperature and a uniform density of 
the nuclei may improve this agreement. 

Similar results are obtained when we apply the model to the Raa measured [23] for 
the D-mesons at the LHC (Fig. [2]). We again see that the model using Ap oc p provides a 
good description of the data over the entire pt range, while that using Ap oc Jp seems to 
describe the data for pt > 4 GeV/c. The constant momentum transfer collision misses 
the shape of the pt distribution completely though it is able to bracket the numerical 
values over a very narrow p? range of 3-5 GeV/c (Fig. [2]). The predictions for t> 2 axe 
given for a ready reference. 

Several factors could affect the value of the energy loss coefficient 'a' for a given 
mechanism. We have verified that increasing (decreasing) the mean free path, A, by 
0.5 fm results in a decrease (increase) of the coefficient, as such that a_e/A remains 
unaltered. 

We have kept a s fixed at 0.3 while estimating the initial charm distribution. Taking 



the renormalization scale as C ' \jp\ + Mq, with C— 1 or 2 leads to a decrease in the 
value of Raa by 7-10 %, which can then be offset by decreasing a B by about 12 %. 
It is interesting to recall that a Focker-Planck equation given by 



01 = A 

dt dpi 



Mp)f + — (%(p)/) 



(14) 



d_ 

d Pj 

describes the evolution of the distribution '/\ °f charm quarks propagating in quark 
gluon plasma [25] and losing energy due to multiple soft scatterings with light quarks 
and gluons. This leads to a drag, Ai(p), and a diffusion -Bjj(p) on the momentum of 
the charm quark. Assuming that Ai(p) depends on momentum only, we have 

A t (p) = A(p 2 ) Pl , (15) 

and the energy loss dE/dx can be related to drag coefficient A(p 2 ) by 

^ = -A(p 2 )p, (16) 
dx 

where E is the energy of the charm quark, and p its momentum. Considering the average 
temperature of the plasma attained at RHIC as ~ 220 MeV [26] . we can read the drag 
coefficient from Fig. 3(a) of Ref. [27J as ~ 0.02 fm -1 for pt up to 5 GeV/c. We can 
re-write the above equation as 
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Comparing this with one of the ansatzes used for energy loss per collision in the present 
work, namely Ap = ag p, we can write that 
dp a B 

T x = ~- p ' (18) 

Thus the effective drag coefficient , A eS (p 2 )\ can be written as 

A eS ( P 2 ) = ^ | (19) 

which reduces to ct^/A for large values of p 

We thus note that the effective drag at RHIC energies is about 0.04 fm -1 compared 
to 0.02 fm -1 estimated for soft multiple collisions by authors of Ref. [27]. Thus we 
conclude that at RHIC energies only half of the energy loss could be due to collisions, 
while the other half could attributed to radiations of gluons. Similarly estimating the 
average temperature at LHC as about 270 MeV and using the results for drag due to 
collisions as ^0.04 fm -1 from Ref. [27] at high momentum, we note that the collisions 
account for only one-third of energy loss at 2.76 TeV/nucleon. 

It is of interest to compare our results with other studies on medium modification 
of charm propagation reported in the literature. Thus Moore and Teaney [28] have 
calculated the diffusion 'D' and drag coefficient 'r^', (denoted by 'A here) using LO 
pQCD as 

D»^(°4)' (20) 



2ttT V 
and 

" D = 4s (21) 

Taking a s ~0.3, this provides T]d ~0.06 fm -1 at 220 MeV and rjo ~0.09 fm -1 at 
270 MeV, which are larger than our values at RHIC and smaller than those at 2.76 
TeV/nucleon at LHC. 

Recall however, the results of Bass et al. [22] that a description of medium 
modification as well as v% for single electrons at RHIC brackets the diffusion coefficient 
'£)' between 1.5/2ttT and 6.0/27rT, when flow and contributions of bottom electrons 
is included. This provides a large value for the drag coefficient between 0.17 and 0.68 
fm -1 . At first these large values may look surprising. However from Eq. [20] we see that 
these would correspond to values of a s between psI.O and ~0.5, which are rather large 
and make the use of perturbative QCD questionable. 

In this connection the studies of Gossiaux et al [30], are also of considerable 
interest which suggest that the collisional energy loss could be substantially larger if 
the Debye mass is replaced by hard thermal loop calculation and a running coupling 
constant is used. We further recall the work of authors of Ref. [31], which suggests 
that a considerable drag could be produced by the resonant heavy quark-light quark 
interaction, beyond that determined by LO pQCD interactions. 

We have already mentioned that c quarks will materialize mostly as D-mesons. 
This can have an interesting consequence which is already apparent in the Figs. [1] and 
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[2J The c quarks after losing energy will pile up at lower energies and this would result 
in a characteristic increase in the production of .D-mesons as well as single electrons 
having lower transverse momenta (Fig. [3]). We note that while there is suppression by 
almost a factor of 2-4 depending on the incident energy (for single electrons) at larger 
Pt, there is essentially no suppression at lower p? at RHIC and even an increase by a 
few percent at LHC where the energy loss is higher and the momentum spectra of c 
quarks have less steeper slopes. The increase in the case of .D-mesons at LHC is rather 
spectacular. This is different from the normal enhancement of mesons having low pt due 
to the so-called Cronin effect, which is expected to be less important at higher incident 
energies and heavy mesons (see e.g., Ref. [32]). We recall that a similar enhancement 
in D-meson production at low p? has been predicted by considering the drag suffered 
the heavy quarks in the plasma and during their hadronization [33J. In a forthcoming 
paper we shall show that this can lead to a slight increase in the production of low mass 
dileptons due to correlated charm decay [34] . 

At the very minimum the present work describes a simple procedure to implement 
energy loss of heavy quarks in relativistic collision of heavy nuclei. It will be of interest 
to explore the energy and centrality dependence of the momentum loss coefficient a. 

4. Summary 

We have used a simple model where charm quarks traversing QGP lose energy in 
multiple collisions. Exploring different parametric forms of energy loss we find that 
the form where the momentum loss per collision is proportional to the momentum 
gives a good description of nuclear modification of single electron production at RHIC 
and D-meson production at LHC. Comparing our results with drag coefficients due to 
collisional energy loss calculated earlier, we find that only a part of the energy loss could 
be due to collisions. We note a characteristic increase in production of single electrons 
and D-mesons having low momenta due to energy loss of charms having larger energy. 
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